Dephasing of a superconducting flux qubit 
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In order to gain a better understanding of the origin of decoherence in superconducting flux 
qubits, we have measured the magnetic field dependence of the characteristic energy relaxation time 
(Ti) and echo phase relaxation time (r|'^'^°) near the optimal operating point of a flux qubit. We 
have measured TJ'^''" by means of the phase cycling method. At the optimal point, we found the 
relation T|^''° 2Ti. This means that the echo decay time is limited by the energy relaxation (Ti 
process) . Moving away from the optimal point, we observe a linear increase of the phase relaxation 



rate (1/T|'^''°) with the applied external magnetic flux. This behavior can be well explained by the 
influence of magnetic flux noise with a 1// spectrum on the qubit. 



As a consequence of the progress in the micro- 
fabrication technology, in the last years solid-state based 
artificial quantum two level systems (TLS), which can 
be fabricated, operated and read out in a controlled way 
have been successfully realized 0,13 • These artificial TLS 
can be used as quantum bits (qubits), i.e. the fundamen- 
tal devices for quantum information systems |^ . 

Among the solid state based qubits one of the most 
promising candidates are the Josephson junction based 
qubits 111, 0, H . Because superconductivity is a macro- 
scopic quantum phenomenon, superconducting qubits 
have an inherent advantage over qubits based on micro- 
scopic systems such as single atoms 0. Furthermore, 
Josephson junction based qubits also represent promis- 
ing candidates for the study of fundamental aspects of 
macroscopic quantum systems 0. For these studies , 
long coherence time as well as high readout fidelity are 
required. Evidently, the fidelity can be improved by re- 
fining the detection method, such as Josephson bifurca- 
tion amplifier technique 3. However, decoherence is still 
determined by the environment of the qubit. Thus, in or- 
der to reduce decoherence we have to first understand its 
origin in detail. 

In this letter, we present a systematic study of the de- 
phasing in a superconducting flux qubit whose level split- 
ting is sensitive to the applied external magnetic flux. 
This allows us to investigate the dependence of the co- 
herence time on the applied flux bias. It is expected 
that such a measurement gives information on the in- 
terdependence between the qubit decoherence and the 
magnetic fluctuations. We present experimental results 
on the magnetic flux dependence of both the energy re- 
laxation time (Ti) and the echo phase decay time (Tl'^^") 
of a superconducting flux qubit. Based on our data, we 
discuss possible mechanisms responsible for decoherence 
in superconducting flux qubits. As a result, we find that 
at the optimal point the coherence time is limited by the 
energy relaxation (Ti process). Away from the optimal 
point, the phase relaxation time is determined by 1// 
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FIG. 1: (color online), (a) Scanning electron micrograph of 
the sample . Three Josephson junction flux qubit is sur- 
rounded with a SQUID detector loop. We apply MW and 
adiabatic dc shift pulse via the MW-line. The mutual induc- 
tance between the MW-line and the flux qubit is ~ 1 pH. The 
readout pulse is applied via the current bias line and qubit 
state is detected by measuring the voltage through the read- 
out line. We use copper powder filter as low-pass filter (LPF) 
of these lines, (b) Pulse sequence for the Ti measurement. We 
generate composite pulse by adding adiabatic DC pulse and 
MW pulse. The dc-SQUID readout pulse is programmed to 
reach the sample just after the adiabatic shift pulse is turned 
off. A rising time of each MW pulse is faster than 0.5 ns. The 
TV pulse width is about ~2 ns. (c) Pulse sequence for the T|'^''° 
measurement. In order to avoid a reduction of the visibility 
by the Ti process, we programmed the end of the adiabatic 
shift pulse and the start of the readout pulse to be reached to 
the sample as soon as the second n/2 pulse finished. 



magnetic fluctuations. 

The sample consists of a three junction flux qubit 
which is surrounded by a dc SQUID Q . Unlike in other 
qubit designs Q, the SQUID loop does not share a cur- 
rent line with the qubit loop. Thus the coupling be- 
tween qubit and readout SQUID is purely inductive (see 
Fig.nja)) with the mutual inductance ^ 13 pH. The elec- 
tromagnetic environment of the qubit is engineered by 
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equipping the SQUID detector with an on-chip shunting 
capacitor. For the fabrication of the AI-AI2O3-AI Joseph- 
son junctions we used usual angled shadow evaporation 
technique 0- The sample was mounted in a dilution 
refrigerator and cooled it down below 50mK. 

At the optimal flux bias point, the mean value of the 
persistent current vanishes for both states of the qubit. 
To be able to discriminate between the two states of the 
qubit at this operation point, we used the adiabatic shift 
pulse method as illustrated in Fig.^b,c) Then we 
apply an adiabatic magnetic dc-pulse with a rise time of 
0.8 ns to the qubit through the coplanar RF-linc. The 
qubit operation is performed during this adiabatic shift 
pulse. At the end of the shift pulse, the qubit will move 
back to the initial bias point adiabatically, where the 
readout can be done. In this way we can get information 
about the qubit state even close to the optimal point. 

We first performed MW spectroscopy of the qubit. 
The resonant frequency (corresponding to the energy 
splitting of the qubit) is determined by saturating the 
state of the qubit with a long enough MW pulse and 
measuring the SQUID switching probability afterwards. 
Figure. El shows the result obtained for the magnetic 
flux dependence of the level splitting of the qubit. In 
the two-level approximation the qubit Hamiltonian is 
Hqb = \ (e($z;)crz + Aa^;), where e{<^x) is the level split- 
ting induced by the external flux $3; threading the qubit 
loop and ax^ Oz are the Pauli matrices. From the spec- 
troscopy experiment we determined the qubit gap fre- 
quency (i.e. the energy splitting at the optimal point) to 
be A//i = 3.9 GHz, where h is the Planck constant. We 
also found the persistent current away from the optimal 
point, 1^ = /p = 370 nA. 

Next, we measured the magnetic flux dependence of 
Ti and Tf^°. Ti is the energy relaxation time describing 
the transition of the qubit from the excited state to the 
ground state due to the interaction with the environment 
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FIG. 2: (color online) . Magnetic flux dependence of frequency 
sweep spectrum. The qubit spectrum is observed as a hyper- 
bola. The field insensitive resonance at 6.2 GHz is attributed 
to the on-chip LC resonator [lol . The red dotted line repre- 
sents a fit of the qubit spectrum by a hyperbolic function. 
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FIG. 3: (color online). Magnetic flux dependence of Ti (open 
squares) and T^'^^" (fllled circles). Although there is consid- 
erable scatter in the Ti data, obviously there is only a weak 
flux dependence. The largest 7^2'^'^° values are obtained at the 
optimal point. Since the data have been obtained via 

the phase cycling technique, the difference of the magnitude 
of error bars of Ti and ri"^*^" is due to differences in the av- 
eraging number (each value is obtained from roughly 
eight times as many data as the Ti value) 



(for k^T ^ A). In order to measure Ti, we use the pulse 
sequence shown in Fig. H^b). Repeating this sequence 
provides the relaxation probability as a function of the 
waiting time. Ti is obtained by an exponential fit of 
this curve. The 7'°'=^° time is the decay time in the echo 
experiments depicted in Fig.^c). 

To understand the effect of non-ideal control in the 
echo sequence we calculate the state of the qubit after 
applied three pulse sequence in the absence of fluctua- 
tions. The integrals (rotation angles) of the three pulses 
are 61 , 62 and ^3 . The time interval between the first and 
the second pulses is ri, while between the second and the 
third pulse it is T2 . We assume that the initial qubit state 
is along the z-axis of the Bloch sphere and calculate the 
z-axis projection of the final state (dz): 

((T2 ) = cos 61 cos 02 cos 63 — sin 61 sin 62 cos 6*3 cos Sluti 

— cos 9i sin 6*2 cos 63 cos 6ujt2 

— sm9i (cos 6*2 + 1) sin 03 cos {Suj (ti + T2))/2 

— sin 6*1 (cos 6*2 — 1) sin 63 cos {Suj (ti — T2))/%X) 

Here, Sto is the detuning. For an ideal echo sequence 
we have Oi ^ 63 ~ 7r/2, 62 = tt, and ti — T2. Under 
these conditions only the fifth term of Eq. (Q, the so- 
called echo term, survives and it is independent of t = 
Ti +T2j (ti = T2) 0. By using phase modulated pulses 
we can control the rotation axis of the Bloch vector 
for all three pulses of the echo sequence The phase 
cycling method makes it possible to extract only the echo 
term. We are able to obtain stable echo decay curve 
during long time experiment. 
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Figure. Ulshows the magnetic flux dependence of Ti and 
Tl'^'^" in the vicinity of the optimal point. The Ti values 
show some scatter but the flux dependence of Ti is weak. 
On the other hand, the T2^^° values show a pronounced 
flux dependence. In particular, there is a sharp maximum 
at the optimal point. This tendency was also reported for 
the edge shared SQUID-flux qubit by Bertet et al. 
At the optimal point, we have T^'^'^" > Ti. 

It is well known that the energy relaxation con- 
tributes to the dephasing process via 



iccho^ 



= (2Ti)^ 



(2) 



Here is the pure dephasing rate due to the fluctuations 
in the energy splitting of the qubit. The dephasing rate 
F;^ consists of two contributions Ftp — ' ^* 



F" +F* FO should 
be identifled with the non-magnetic dephasing processes 
or with the second-order contribution jlfl] . F* is the first- 
order contribution due to magnetic dephasing processes. 
At the optimal point, the magnetic flux derivative of the 
energy dispersion is flat | | — 0, and there is no con- 
tribution of first-order magnetic flux fluctuations to the 
dephasing rate, i.e. F* = 0. We obtain T°'''^° = 250 ns 
at the optimal point and the average value of Ti around 
the optimal point is about 140 ns. From these values, we 
estimate the pure dephasing rate at the optimal point to 
be F^ = (2.3/is)~^. Remarkably, the dephasing time of 
the qubit (r|'^'^°) is almost entirely determined by the en- 
ergy relaxation process, i. e. T^"^^" « 2Ti. In other words, 
the dephasing of the flux qubit near the optimal point is 
predominantly determined by the high-frequency noise 

For dephasing dominated by magnetic flux noise with 
a smooth spectrum near lu = this rate is given by 

2 



dAE 



9$. 



S.s>{u; = 0) 



(3) 



Here AE = AE{^x) is the external, flux dependent qubit 
energy splitting and S^{uj = 0) is the low frequency 
component of a correlation function of flux fluctuations. 
From AE = ^e2($^) + A2 we obtain dAE/d^^ = 
COST] {de/d^x), where tanrj = A/e. Since the variation 
of the bias with external flux, de/d^x, is usually almost 
constant over a wide interval of the applied flux. F* be- 
comes proportional to cos^ rj. Taking into account that, 
according to our experimental findings, the Ti time is al- 
most independent of the applied flux and that a relation 
cosTy cx e is satisfled near the optimal point, a parabolic 
behavior of l/T^'^^" is expected in this region. However, 
as shown in the inset of Fig. |31 the measured 1/T|'^'^° 
versus applied flux curve is linear. This result is very dif- 
ferent from the prediction of Eq. ((^J. As a consequence, 
our experimental data cannot be explained with a simple 
Bloch-Redfield type of decoherence theory where a noise 
with short correlation time (white noise near w = 0) is as- 
sumed. In order to explain the observed behavior we need 
to consider decoherence from a long correlated noise. 



Even when the pure dephasing is non-exponential the 
echo decay curve can be represented as the product of en- 
ergy decay (Ti) and pure phase decay. In this case Eq. (0 
is replaced by a non-exponential decay curve for the ex- 
pectation value of (Tj at the end of the echo sequence 

p(t) = {a^) = exp (-^It ) • /°ccho [t) , where pccho {t) = 
(exp (§) cosr; $ (r) dr - /| $ (r) dr] }). As- 
suming Gaussian fiux fluctuations this gives pccho {t) ~ 



exp 



In the case of 



white noise, S^{uj) = const., Pccho (0 gives a simple 
exponential decay curve, and the flux dependence of 
the decay rate follows Eq. On the other hand, 

for S^{u) ^ ^ ^ ^ (l//-noise |l3) Pccho (t) = 



exp 



A 
f 



'echo 

This decay law is 



Gaussian and the relaxation rate is proportional to 

We do not attempt to determine the echo decay shape 
from the observed relaxation curve because our experi- 
mental data do not have precision enough to distinguish 
between the different types of decay curves (e.g. simple 
exponential, Gaussian or algebraic). Instead, we assume 
a special decay shape based on the 1// fluctuations spec- 
trum p {t) = exp (-2I7) • exp (-r°t) • exp (^-F* 

And we flnd the magnetic flux dependence of F* by flt- 
ting of this formula. In this fltting Ti and F° are fixed 
and F* is the fitting parameter. This relaxation func- 
tion is constructed from three components. The first one 
describes the energy relaxation (Ti process). In the sec- 
ond component F^ is the phase relaxation rate, which 
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FIG. 4: (color online). Magnetic flux dependence of val- 
ues. The values are derived from fits of the data by the 1// 
fluctuation model using a constant Ti value. The filled circles 
represent F^^ values estimated using an average Ti value of 
140 ns. Since the T\ values have some scattering, we also es- 
timated values using both Ti = 140 ± 30 ns. The triangles 
and inverted triangles represent the values estimated from 
Ti = 110 ns and 170 ns, respectively. The red solid line shows 
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is caused by non-magnetic fluctuations [ig] and second- 
order contributions of magnetic fluctuations. Near the 
optimal point, the second-order contributions of 
is small then we can neglect it. Finally, F* is the phase 
relaxation rate due to 1/f magnetic fluctuations. We ob- 
tain Ti = 140 ns from the energy relaxation curve and 
F° — (2.3 ns)~^ from the echo decay rate at the optimal 
point. Figure. 01 shows the magnetic flux dependence of 
the resulting F*. Very close to the optimal point, the 
contribution of F* to the echo decay curve is very small 
compared to that of 1/Ti. In this region the F* val- 
ues cannot be extracted with sufficient accuracy from 
the echo decay curve. Apart from this region, the mag- 
netic flux dependence of F* is well fitted by expression 
F* oc I I — I M cos 77] . This means that F* is propor- 
tional to absolute value of the magnetic flux difference 
from the optimal point. Thus our experimental result 
suggests the existence of 1// type magnetic fluctuations 
in the frequency range of cj « O (F*) . The magnitude of 
the 1/f flux noise can be estimated for the A value of the 
correlation function 5*$ (uj) — From the spectrum 

dispersion and the magnetic flux dependence of F* we 

obtain A « (lO~^<i>o)^. This value is of the same order 
as for other superconducting circuit Il9|. 

In summary, we have performed a systematic study 
of the phase and energy relaxation of a superconducting 
flux qubit by measuring the magnetic fiux dependence 
of the energy relaxation (Ti) time and of the echo de- 
cay time (T|'^'^°). At the optimal point we find that the 
contribution of magnetic fluctuations is negligible. The 
jiccho i^jjjjg Qf ^j^g qubit is then determined by the energy 
relaxation resulting in J'^'^^o ^ 2X1. Thus, in order to 
improve the coherence time of our flux qubit at the opti- 
mal point, we need to suppress high frequency transversal 
fluctuations Sx^y (A/h). Away from the optimal point the 
magnetic flux fluctuations dominate the echo dephasing 
time Tl^^ho^ rp]-^g observed linear dependence of F* on 
magnetic flux around the optimal point suggests that the 
magnetic fluctuations have a 1/f -type spectrum in a fre- 
quency range of the order of 1/T2'^^°. To our knowledge 
these results provide the first experimental clue for the 
possible origin of decoherence in a superconducting fiux 
qubit, which is not sharing an edge with its dc SQUID 
detector. 
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Note added: After complete this work, we became 
aware of similar results by F. Yoshihara et al. [2(1 |. 
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